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ABSTRACT 

We  consider  a  finite  one-dimensional  elastic  body  with  stored-energy 

2  2 

function  of  the  form  W(Y)  +  e  (y*)  ,  where  Y  is  the  strain,  y 1  the 
strain  gradient,  and  c  a  small  parameter.  We  assume  that  W(Y)  is  non- 
convex,  of  the  type  capable  of  supporting  two-phases.  We  show  that  when  the 
body  is  acted  on  by  end  loads  the  only  stable  strain  fields  are  the  single¬ 
phase  solutions  Y  =  constant. 
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SIGNIFICANCE  AND  EXPLANATION 


Cf  '  ,  't  r,  ?uk  2 


^  n  i~  S.^,i'kar< 

-Wtf  consider  a  one-dimensional  elastic  bar  acted  on  by  end  loads.  To 
allow  for  the  possibility  of  phase  transitions  we  consider  a  stored-energy 
function  of  the  form  W(y)  +  (f /(Y*)  with  W(^f >  a  nonconvex  function  of 

<L 

strain  jf,  a  function  of  the  type  capable  of  supporting  two-phases.  The  term 

^  1  '*•  N*  i. 

involving  the  strain  gradient  is  added  to  model  the  region  of  rapidly 

varying  strain  in  the  vicinity  of  a  phase  transition.  „  . 

_  l:  p  f  •  !o  « 

We  showr^that  this  problem  has  smooth  solutions  which,  for  t  small, 

correspond  to  the  single-phase  and  two-phase  solutions  of  the  problem  with 

C „  <■ 

f  ■  0.  We  shown^urther  that  the  only  stable  solutions  of  this  problem  are  the 
single-phase  solutions  Y  =  constant. 


! 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
susnary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ONE- DIMENSIONAL  STRUCTURED  PHASE  TRANSFORMATIONS 
UNDER  PRESCRIBED  LOADS 


Jack  Carr,  Morton  E.  Gurtin,  and  Marshall  Slemrod 
1.  Introduction 

Consider  a  bounded  one-dimensional  body  identified  with  the  interval 

2 

[ 0,L] .  We  assume  that  the  body  is  elastic  with  stored  energy  W(Y)  a  C 
function  of  strain*  Y»  0  <Y  <  *,  end  with  stress 


o'  >  0  on  (0,y  ) ,  o'  <  0 

* 

on  (y  ,Y,)»  end  o'  >  0 

on  (Y *»“)•  Let 

*  * 

-  o(Y*),  o  *  o (Y  ). 

# 

Then  for  f  e  [ot,a  ) 
there  are  exactly  three 
values  of  y  for  which 

o(y)  *  f  • 

We  denote  these  by 

a(f )  <  1(f)  <  0(f)  , 

omitting  the  argument  f  when  convenient,  and  we  extend  a(f)  {respectively 

* 

0(f)}  in  the  obvious  manner  for  f  €  (0,o4J  {respectively  [o  ,•»)}. 


o(Y) 


of  the  form  shown  in  Figure  1.  Thus 


dW(Y) 

dY 


Figure  1.  The  stress-strain  law. 


1 

If  the  underlying  deformation  u  carries  material  points  x  e  (0,L)  to 
points  u(x),  then  y  ”  du/dx. 
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We  restrict  our  attention  to  loading  applied  by  a  soft  device  ,  so  that 
the  load  f  («  constant)  is  prescribed.  The  problem  of  determining  the  stable 
strain  field  resulting  from  this  load  consists  in  minimizing  the  corresponding 
potential  energy.  This  problem  can  be  stated  precisely  as  follows: 

(5)  Minimize 

fg  [w m  -  fyldx  (1.1) 

over  all  Y  e  L^O.L)  with  W(Y)  e  lVo.L). 

This  problem  has  been  studied  by  Ericksen  [1975],  and  we  briefly  record 
some  of  his  results.  The  Euler-La grange  equation  is 


Figure  2.  The  Maxwell  line. 


o(Y)  *  f  / 

while  the  form  of  the 
minimizer  depends  on 
whether  or  not  the 
load  f  is  smaller  than, 
equal  to,  or  larger  than 
the  Maxwell  stress 
oM»  where  0M  is  the 
unique  value  of  stress 
for  which  (Figure  2) 


The  strains  a  and 
H 

respectively. 


W(B  )  -  W(a  )  »  a  (B  -a  ) 

rl  Pi  H  PI  PI 


°M  *  a(<V'  SM  “  B(V 


Bm  are  the  lower  and  upper  Maxwell  strains , 


2 

Cf.  Ericksen  [1975]. 


I 
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Figure  3.  The  potential  W(y)  -  fy  • 


As  Is  clear  from  Figure  1,  the  potential  W(y)  -  fy  has  the  fora  shown 

in  Figure  3.  For  f  <  the  local  alniaua  at  8  disappears,  as  does  the 

* 

local  miniaua  at  a  when  f  >  o  .  Therefore  for  f  f*  the  solution  to 
Problem  S  is1 

y(x)  =  a  when  f  <  <»M  , 
y(x)  =  8  when  f  >  aM  • 

for  f  «  c  the  problea  is  far  aore  interesting.  Here  any  function  of  the 

H 

fora 

Y(x)  ■  a  for  all  x  6  S  , 
n 

y(x)  -  8m  for  all  x  6  tO,L]\S  ,  (1.2) 

S  c  [0,L] ,  S  measurable  , 

is  a  ainiaiser. 

We  will  use  the  following  terminology »  solutions  of  the  fora 

y(x)  =  constant 

will  be  called  single-phage  solutions >  solutions  of  the  fora  (1.2)  with  S 

_ 

The  energy  (1.1)  also  has  a  local  miniaua  relative  to  L  (0,L)  at  a  for 
<  f  <  a*  and  one  at  8  1 for  a*  <  f  <  omj  these  constant  functions  are  not 
local  minima  relative  to  L  (0,L). 


-3- 


and  [0,^X8  nontrivial  will  ba  called  two-phaaa  aolutlonai  in  the  latter 
case,  when  both  the  nuaber  n^  of  nontrivial  connected  components  of  8  and 
the  number  n2  of  nontrivial  connected  coaponenta  of  [0,10X8  are  finite/ 
we  will  refer  to  n^+n2-1  as  the  degree  of  the  solution. 

Thus  for  f  f  o  the  solution  is  single-phase  (and  unique)}  for 

H 

f  *  a  there  are  two  single-phase  solutions  (Y  =  <*  and  y  =  6  )  and  an 

uncountable  infinity  of  two-phase  solutions  of  each  degree  >  1. 

For  the  case  f  -  o  it  seems  reasonable  to  ask  whether  any  of  the 
n 

solutions  are  preferred.  One  might  expect  that  the  single-phase  solutions 
are,  in  soae  sense,  more  stable  than  the  two-phase  solutions.  We  here  show 
that  this  expectation  is  indeed  true  provided  we  consider  the  present  theory 
as  an  approximation  to  a  higher-order  theory  which  allows  the  strain-gradient 
to  enter  the  constitutive  equation  for  the  stored  energy. 

A  general  nonlinear  theory  of  this  type  was  apparently  first  given  by 
TOupin  [1962]  and  is  based  on  constitutive  functions  for  the  stored  energy  and 
stress  of  the  forms 

A 

W(Y,Y')  and  o( Y,Y’> 

with 

o(Y,Y')  -  W(Y. Y’>  • 

The  presence  of  the  strain-gradient  Y'  *  dy/dx  models  situations  in  which  the 
strain  varies  rapidly;  in  this  instance  the  force  interactions  within  the 
underlying  crystal  lattice  cannot  be  adequately  described  by  the  usual  notion 
of  stress  and  one  needs  -  in  addition  to  the  (classical)  stress  a  -  a  second- 

A 

order  atress  u  which,  in  the  present  theory,  is  related  to  W  through 

P(Y,Y,>  -  W(YfY')  • 
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1 

Energy  functionals  of  this  type  originate  in  the  classic  paper  of  van  der 
ffaals  [1893].  Cahn  and  Hilliard  [1958]  independently  rederived  van  der  Waal s' 
basic  theory,  and  in  the  last  twenty  five  years  gradient  theories  have  bee one 
a  popular  tool  in  analyzing  phase  transitions  and  other  physical  phenomena. 
(Cf.  Rowlinson's  [1979]  translation  of  van  der  Waal  [1893]  for  a  list  of 
selected  references.) 
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theory  to  apply,  and  for  that  reason  the  boundary  conditions  (1.5)  seen 
appropriate.  [Of  course,  one  might  expect  it  to  be  difficult  if  not 
impossible  to  achieve  this  form  of  boundary  condition  in  practice.] 

It  is  clear  that  if  a  single-phase  solution  Y  minimises  the  original 
functional  (1.1),  Y  trivially  minimises  (1.3),  as  Y'=  0.  Thus  the  complete 
list  of  solutions  of  Problem  S£  consists  of  the  single- phase  solutions 


Y(x)  S  o 


Y(x)  =  oM 

yu)  =  eM 


} 


for  f  <  , 

for  f  -  oM  , 


Y(x)  s  6  for  f  >  oM  . 

2  2 

Because  of  the  term  e  ( y " )  in  (1.3),  two-phase  solutions  of  the  form 

(1.2)  are  not  members  of  any  function  class  which  renders  S a  meaningful 

problem.  On  the  other  hand,  there  are  solutions  of  the  Euler-Lagrange 

equation  (1.4)  for  f  »  o  which  spend  most  of  their  time  near  a  and  6 

h  n  H 

and  which  correspond  to  the  two-phase  solutions  of  the  unstructured  theory. 
This  raises  the  question  of  whether  such  structured  solutions  can  be  relative 
minimizers  of  (1.3).  He  will  show  that  the  answer  to  this  question  is  no;  in 
fact,  we  will  prove  that  all  nonconstant  solutions  of  the  problem 

2e2Y"  -  o(Y)  -  f  , 

(1.6) 


Y'(0)  -  Y’(b)  -  0 


are  unstable . 

Remark.  Our  analysis  is  also  appropriate  to  a  fluid  in  a  one-dimensional 
chamber  acted  on  by  a  pressure  p  ■  -f.  In  this  instance  y  is  the  specific 
volume,  W(y)  is  the  internal  energy  per  unit  mass,  x  e  [0,L]  labels 
material  points  in  a  fixed  reference  configuration  with  uniform  density  pQ, 
and  the  functionals  (1.1),  (1.3)  -  multiplied  by  p^  -  represent  the  total 
potential  energy. 


i 
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2.  Structured  solutions. 

It  is  convenient  to  change  length  scales.  We  let 

t  -  x/e,  Le  -  h/e 

and  consider  Y(t)  as  a  function  of  t  rather  than  xj  the  problem  (1.6) 
then  takes  the  fora 


with  Y  *  dY/dt, 
"first  integral" 


<v 


2if  -  0(Y>  -  f  , 


Y(0)  “  Y(Le)  -  0 


(2.1) 


and  multiplying  ( 2 . 1 )  by  Y  and  integrating  leads  to  the 


where 


Y2  -  *(Y)  -  «(Y0)  , 

*(Y>  -  W(Y)  “  *Y  . 
Y0  “  Y(0)  . 


(2.2) 


(2.3) 


For  future  use  we  note  also  that  the  associated  energy  (1.3)  is  now  given  by 

E(Y )  -  J0G  tY2  +  *(Y)ldt  .  (2.4) 

Remark.  Equation  (2.1).,  is  the  equation  of  motion  of  an  undamped  oscillator 
of  "mass"  2  and  "potential  energy"  -*(Y)>  the  functional  E  is  the  corre¬ 
sponding  "Hamiltonian". 

We  now  state  two  propositions!  these  allow  us  to  concentrate  on 

*  * 

nonconstant  solutions  too  Problem  S£  and  to  f  €  (o#,o  ).  The  first 

A 

proposition  follows  from  the  fact  that  Y  =  constant  solves  S£  if  and  only 
if  c(Y)  ”  f.  The  second  is  a  direct  consequence  of  the  phase  portraits  for 

i  *  • 

f  f  (o#,o  )»  for  such  f  there  are  no  trajectories  with  Y  ”  0  at  two 
distinct  places. 
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Proposition. 

follows: 


Ths  complete  list  of  constant  solutions  of  Problem  5  is  as 


a,X, 

M. 

B, 


0  <  f  <  o#  , 


o#  <  f  <  a 


(2.5) 


f  >0 


Proposition.  All  nonconstant  solutions  to  Problem  S£  correspond  to  f  in 

* 

the  range  (o#,o  ). 

* 

The  graphs  of  -*  for  f  e  (o#,o  )  are  simply  the  curves  of  Figure  3 
turned  upside-down.  In  view  of  the  remark,  this  allows  us  to  easily  arrive  at 
the  phase  diagrams  shown  in  Figure  4. 


Figure  4.  Phase  diagrams 


Let  us  agree  to  use  the  term  admissible  trajectory  for  a  nonconstant 
solution  Y(t),  0  <  t  <  T^,  of  the  problem 

2y  -  o(Y)  -  «  # 


r 


e  (0,*)  is  then  the  duration  of  Y,  and  the  number  of  times  t,  not 


including  t  ■  0,  that  Y(t)  *  0  is  the  degree  of  Y*  Problem  S£  (for 
nonconstant  solutions)  can  then  be  stated  as  follows:  find  an  admissible 
trajectory  with  duration  L£. 

this  formulation  allows  us  to  easily  establish  the  following  existence 

A 

theorem  for  Problem  Se  • 

* 

Theorem .  Let  f  S  (ot,o  ).  Then  there  is  a  constant  C  »  C(f)  >  0  such  that 
(i)  for  e  >  L/C,  Problem  Se  has  only  constant  solutions ; 

A 

(ii)  for  e  <  L/C,  Problem  has  nonconatant  solutions  of  degree  n  for 
each  n  ■  1,2, . . . ,N ,  where  H  i£  the  largest  integer  with  N  <  L/( ce ) . 
Proof.  We  will  establish  this  result  only  for  the  case 


the  remaining  cases  o  <  f  <  a  and  a  <  f  <  a  are  completely  analogous. 

It  is  clear  from  the  phase  portrait  that  each  admissible  trajectory  Y 

may  be  labeled  by  its  initial-value  Y(0)  *  YQ  and  its  degree  N.  Further, 

each  such  Y  has  y  in  either  (a  ,X  )  or  (X  , B  ),  and  there  is  an 

U  MM  MM 

obvious  one-to-one  correspondence  between  the  two  resulting  classes  of 
trajectories;  this  allows  us  to  restrict  our  attention  to  trajectories  Y 
with 


Given  such  a  trajectory,  Y(t)  >  0  for  0  <  t  <  tj,  where  t1  =“  T^/tl  with 
the  duration  and  N  the  degree  of  yt  further,  Y(t1>  “  Y^  *  where 

is  the  unique  solution  of 

*(Y1 >  -  *(Y0) 

in  the  interval  (*„#&„)  v’cf.  (2.3)  and  Figure  3).  The  duration 
M  M 

Ty  "  T<VN> 
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of  such  a  trajectory  is  then  trivially  equal  to 

fy*  •  -1 

Ty  -  *  /0T  y(t)  <mt)  , 

and  if  we  change  variable  of  integration  from  t  to  y»  we  conclude,  with 


the  aid  of  (2.3),  that 


T(Y0,H)  -  WP(Y0)  , 


fW 

T(Yfl)  “  /  ' _ 

T0  /*(C)-*(Y„) 


A L 


(2.6) 


T(Yq)  has  the  limits 


T(a 


•  ♦  -•  *V  -  v5frr > 0 


(2.7) 


(the  latter  being  one-half  the  period  of  the  oscillator  defined  by  (2.1)^ 
linearized  about  the  rest  point  X^) •  Further  the  integrand  in  (2.6)  is 
bounded  away  frcm  zero,  as  is  the  interval  of  integration  for 
Yg  8  Thus*  in  of  (2. 7 >2»  C  defined  by 


C  -  inf{T(Y0):  ^  <  Y„  < 


is  strictly  positive. 

Clearly,  the  equation 


T(Yq>  -  L£  (2.8) 

has  no  solution  Y-  S  (a ,A  )  when  L  <  C,  and  this  implies  (i)  of  the 

U  MM  t 

theorem. 

On  the  other  hand,  for  L£  >  C  we  conclude  from  (2.7) ^  that  (2.8)  has  at 

least  one  solution  Yn  8  (o  ,X  ),  as  does  the  more  general  relation 

U  MM 

nT(Yfl)  "  Le 

for  n  »  1,2, . . . ,N  with  M  the  largest  integer  <  L£/C.  This  yields  (ii). 


1 

Cf.,  e.g..  Hale  [1969],  p.  179;  Arnold  [1978],  Chapter  2,  Section  4. 
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Remark.  Let  f  ■  o  and  consider  the  solution  of  degree  1  starting  in 

<«  ,X  ).  For  e  small  this  solution  will  be  close  to  the  heteroclinic  orbit 

from  a  to  0  and  hence  will  spend  most  of  its  time  near  a  and  B  •  in 
n  n  mm 

this  sense  the  solution  represents  the  structured  analog  of  the  unstructured 
(€*0)  two-phase  solutions  of  degree  1  which  have  a  followed  by  B  . 

M  M 

Similarly,  the  structured  solution  of  degree  2  starting  in  (a  .A  )  iB  the 

11  M 

analog  of  the  unstructured  two-phase  solutions  of  degree  2  which  have  a 

M 

followed  by  B^  followed  by  a^.  Similar  assertions  apply  to  the  other 
nonconstant  solutions. 


i 


3.  Stability. 

Throughout  this  section  the  energy 

E<Y>  -  /0E  {Y2  +  W(Y>  -  fYldt 
is  considered  as  a  functional  with  domain 

r  -  (Y  e  H^O.L  ):  Y  >  0  on  [0,1] }  . 

In  addition,  we  assume  that  W  is  a  C  function,1  so  that,  by  (2.1)^, 

a 

solutions  of  Se  are  C  on  [0,Le]. 

Let  Y  be  such  a  solution.  Then 
( i )  Y  is  globally  stable  if 

E(Y>  <  E(<d)  for  all  u  e  T  » 

(ii)  Y  is  locally  stable  if  there  is  a  neighborhood  fl  of  Y  in  T  such 
that 

E(Y )  <  E<(i>)  for  all  w  6  (1  i 

OD 

(iii)  Y  is  unstable  if  given  any  norm  on  the  space  C  [0,L£]  and  any 
neighborhood  0  of  Y  in  this  normed  space  there  is  an  u  e  Q  with  lower 
energy : 

E(u)  <  E(Y>  . 

2  * 

Theorem.  Ml  nonconatant  solutions  of  Problem  Sc  are  unstable . 

The  proof  is  based  on  the  following  leasu,  which  is  well  known. 

*  m 

La— Let  Y  be  a  solution  of  $  •  Suppose  there  exist*  en  nec  [0#L  ] 
such  that 

L 

KY.n)  -  J0e[2n2  +  w"<Y)n2]dt  <  o  . 

Then  Y  if,  unstable. 


C2  would  suffice  provided  c'tO^  ]  in  (iii)  is  replaced  by  C2[0,L  ]. 

2  e  e 

Within  a  slightly  different  context  a  similar  theorem  has  been  established  by 
Chafee  (1975].  <Cf.  Casten  and  Holland  [1978]  and  Matano  [1979].)  We  arrived 
at  our  results  before  discovering  the  above  references. 
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V 


Proof.  Lot 


*  (t)  -  e(y  +  tn) 

for  all  sufficiently  small  t.  Then  a  simple  calculation,  based  on  (2.1), 
shows  that 

(0)  -  o,  *"<0>  -  i(Y,n)  . 

Since  I(Y,h)  <  0,  this  tells  us  that  <fi  has  a  strict  relative  maximum  at 
t  •  0/  hence 

E(Yt)  <  E(Y),  Yt  -  Y  +  Th 
for  all  sufficiently  small  t  )<  0.  Moreover,  as 

*yt  -  yi  -  I t |  ini 

(!•!  arbitrary),  every  1*1  -  neighborhood  contains  functions  Y  ,  T  /  0, 
and  the  proof  is  complete. 

A 

Proof  of  Theorem.  Let  Y  be  a  nonconstant  solution  of  S  ,  and  let 

n(t)  -  Y(t)  +  6$(t) 


with 


and  6  an  arbitrary  constant.  Then,  writing  Ig 

2 

Ig  ■  A  +  26b  +  6  c  , 


where 

A  -  J0Et2Y2  +  W"(Y)Y2]dt 


B  -  Joe[2Y§  +  W"(Y)Y31dt 
C  -  I(Y,3)  . 


I(Y,h), 


(3.1) 


(3.2) 


and  hence,  using  (2.1) 


A  -  2  4C(V  Y)*dt  -  0  , 

L 

B  -  2  /Q  (Y  8)  *dt  *  -2Y(0) 


'3.3) 


Clearly,  Y(0>  ¥  0,  for  otherwise,  by  (2.1),  y  would  be  constant)  thus 
B  ¥  0.  We  therefore  conclude  from  (3.2)  and  (3.3)  that  for  some  6^0, 

h  <  0  ' 

and  the  proof  is  complete . 

Remark.  It  is  interesting  to  note  that  the  last  theorem  is  independent  of  the 
particular  form  of  the  stored  energy  W. 

To  complete  our  study  we  have  only  to  investigate  the  stability  of  the 
constant  solutions  (2.5);  but  this  is  not  difficult. 

Proposition.  The  constant  solutions  (2.5)  have  the  following  properties) 

a  is  globally  stable  for  0  <  f  <  a  and  locally  stable  for  o  <  f  <  o  »  8 
Ti_  —  n  "  M 

is  globally  stable  for  o  <  f  <  •  and  locally  stable  for  o#  <  f  <  a; 

* 

X  i£  unstable  for  o#  <  f  <  0  . 

Proof.  Let  Y  be  a  constant  solution.  Then  for  any  u  e  T, 

L 

E(o>)  -  E(Y)  -  /„  (#(«)  -  4<Y)  +  <•»  Jdt 

and  the  global  stability  of  a  follows  from  the  fact  that  •  has  a  global 

* 

minimum  at  a  for  0  <  f  <  a  •  Also,  for  o  <  f  <  o  ,  a  is  a  local  minimum 

M  H 

for  ♦  ;  hence  there  is  a  6  >  0  such  that  4(a)  <  *(t)  and  t  >  0  for  all 
t  with  |t-a|  <  6.  Further 

n  -  {(I)  e  Ti  |u(t)-a|  <  6  for  all  t  C  [0,LeJ} 
is  a  neighborhood  of  Y(t)  =  a  in  h\o,L£)  and  E(a)  <  E(o>)  for  all 
<o  e  fl;  hence  a  is  locally  stable  for  oM  <  f  <  o  . 

The  assertions  concerning  8  are  established  in  a  similar  manner. 
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* 


Finally,  for  o#  <  f  <  o  ,  X  la  not  a  local  minimum  for  ♦  and  we  can 
find  points  t  which  are  arbitrarily  close  to  X  and  have  *(t)  <  *(X). 

m 

Thus  we  can  find  constant  functions  u  which  -  in  any  norm  on  C  [0,1.^]  ~  are 
arbitrarily  close  to  Y(t)  =  X  and  have  E(o>)  <  E(y).  Therefore  X  is 
unstable . 
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